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, ABSTRACT 

The dust feature G159.6-18.5 in the Perseus region has previously been observed with 
^ the COSMOSOMAS experiment (Watson et al. 2005) on angular scales of ^ 1° and 

^ was found to exhibit anomalous microwave emission. We present new observations of 

CS| this dust feature, performed with the Very Small Array (VSA) at 33 GHz, to help 

OO increase the understanding of the nature of this anomalous emission. 

On the angular scales observed with the VSA (^ 10 - 40'), G159.6-18.5 consists 
of five distinct components, each of which have been individually analysed. All five 
^\ of these components are found to exhibit an excess of emission at 33 GHz, and are 

found to be highly correlated with far-infrared emission. We provide evidence that 
^\ each of these compact components have anomalous emission that is consistent with 

electric dipole emission from very small, rapidly rotating dust grains. These compo- 
^ nents contribute ~ 10 % to the fiux density of the diffuse extended emission detected 

•i-H by COSMOSOMAS, and are found to have a similar radio emissivity. 

^ Key words: radiation mechanisms: general - ISM: individual: G159.6-18.5 - clouds - 

dust, extinction - radio continuum: ISM 



1 INTRODUCTION 

Within our own Galaxy, there are three well defined sources 
of diffuse Cosmic Microwave Background (CMB) foreground 
continuum emission, namely, synchrotron and free-free 
emission, which dominate at frequencies below ^ 60 GHz, 
and thermal (vibrational) dust emission, which dominates 
at higher frequencies. However, there is increasing evidence 
supporting the existence of a new continuum emission mech- 
anism in the frequency range ?^ 10 - 60 GHz (Leitch et al. 
1997; Banday et al. 2003; Dickinson et al. 2006; Davies et 
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al. 2006; Hildebrandt et al. 2007; Dobler & Finkbeiner 2008, 
and references therein). 

Various theories have been proposed to explain this 
"anomalous" emission, including magnetic dipole emis- 
sion (Draine & Lazarian 1999), hard, flat spectrum, syn- 
chrotron emission (Bennett et al. 2003) and bremsstrahlung 
from very hot (Te ~ 10^ K), shock heated gas (Leitch et 
al. 1997). However, the currently favoured emission mecha- 
nism is that of very small (N ^ 10^ atoms), rapidly rotating 
(~ 1.5x10^° s~^) dust grains, emitting electric dipole ra- 
diation, commonly referred to as "spinning dust" (Draine 
& Lazarian 1998; Ali-Hai'moud et al. 2009). Fulleranes 
have also been proposed as the source of anomalous emis- 
sion (Iglesias-Groth 2005). 
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Figure 1. ■'^^CO integrated intensity false colour image (Ridge et al. 2006a) of the Perseus molecular cloud overlaid with contours of the 
IRIS 100 iiuv emission (10, 20, 40, 60 and 80 % of the peak emission 680 MJy sr-^) showing the location of the dust shell G159. 6-18.5 
(dashed line), the central star HD 278942 (cross) and the location of 2 major star formation sites, IC 348 and NGC 1333. The box 
illustrates the region in which the VSA observations were performed. 



This anomalous emission has been detected on large 
angular scales across substantial areas of the sky (Kogut et 
al. 1996; de Oliveira-Costa et al. 1997, 1999, 2002; Miville- 
Deschenes et al. 2008) and also on small angular scales in 
specific Galactic objects, such as LDN 1622 (Finkbeiner et 
al. 2002; Finkbeiner 2004; Casassus et al. 2006; Dickinson 
et al. 2007), G159.6-18.5 (Watson et al. 2005), RCW 175 
(Dickinson et al. 2009), LDN 1111 (AMI Consortium; Scaife 
et al. 2009a) and a number of Lynds Dark Nebulae (AMI 
Consortium; Scaife et al. 2009b) and is often found to be 
tightly correlated with the far-infrared (FIR) emission. The 
FIR emission is generally interpreted as due to three dust 
components: big grains (BGs), very small grains (VSGs) and 
polycyclic aromatic hydrocarbons (PAHs), deduced from 
IRAS observations (Desert et al. 1990). 

One of these Galactic objects, G159.6-18.5, is located 
within the Perseus molecular complex. Observations per- 
formed with the COSMOSOMAS experiment (Watson et 
al. 2005) identified a bright source of dust-correlated emis- 
sion at 20 - 30 GHz, with a peaked spectrum, indicative 
of spinning dust. This is perhaps the best example of spin- 
ning dust emission where the peaked spectrum is well fitted 
by a Draine & Lazarian (1998) spinning dust model. These 
results prompted further investigations of this region, which 
were performed with the Very Small Array (VSA), a 14- 
element interferometric array operating at 33 GHz; Watson 
et al. (2003) give the nominal VSA setup. In the present 
paper we present and discuss our results from these new ob- 



Table 1. Summary of the characteristics of the VSA in the super- 
extended array configuration. 



Location 


Teide Observatory 


Altitude 


2400 m 


Latitude 


+28^ 18' 


Declination range 


-7^ < Dec < +63^ 


No. of antennas (baselines) 


14 (91) 


No. of correlations 


182 


Tsys (K) 


^ 35 


Frequency (GHz) 


33 


Bandwidth (GHz) 


1.5 


Primary beam FWHM (arcmin) 


72 


Synthesized beam FWHM (arcmin) 


^ 7 


Sensitivity (mJy beam~^ in 50 hrs) 


^ 5 



servations. Section 2 describes the region under investigation 
and presents the VSA observations, while Section 3 details 
the data reduction. Section 4 discusses the available ancillary 
data, which were used to help gain a better understanding 
of the physical conditions within the region. In Section 5 we 
produce spectral energy distributions (SEDs) for 5 features 
in the region. In Section 6 we analyse these SEDs, investi- 
gate the correlation between the microwave and IR emission, 
and analyse the relationship between our VSA results and 
the COSMOSOMAS results (Watson et al. 2005). Finally, 
our conclusions are discussed in Section 7. 
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2 VSA OBSERVATIONS OF G159.6-18.5 
2.1 The G159.6-18.5 Region 

The Perseus molecular complex is a giant molecular cloud 
located in the Perseus constellation at a distance of 
- 260 pc (Cernicharo, Bachiller & Duvert 1985). The cloud 
chain is ~ 30 pc in length and contains many well known 
sites of active star formation. The dust feature observed 
by Watson et al. (2005), G 159.6-18.5, appears as a remark- 
ably complete shell of enhanced FIR emission with diame- 
ter ^ 1°.5 as shown in Fig. 1. 

Fig. 1 displays a ^^CO integrated intensity image (Ridge 
et al. 2006a) illustrating the extent of the Perseus molec- 
ular cloud. Overlaid on the image are contours of the 
IRIS^ 100 fim emission showing the location of the dust 
shell G159.6-18.5, which does not appear to be traced by 
the ^^CO integrated emission. Also overlaid is a rectangle 
showing the coverage of the VSA observations. 

Due to its shape, G159.6-18.5 was initially believed to 
be a supernova remnant (Pauls & Schwartz 1989; Fiedler 
et al. 1994), but more recent observations have determined 
that the source of the sheh is the O9.5-B0V star, HD 278942, 
at its geometric centre, and that the shell is filled with Hii 
gas (Andersson et al. 2000). Observations performed as part 
of the COMPLETE Survey (Ridge et al. 2006b) suggested 
that G 159. 6-1 8. 5 was indeed an expanding Hll bubble lo- 
cated on the far side of the molecular cloud. We will discuss 
the distribution of the ionized gas in more detail in Sec- 
tion 4.1. 



2.2 VSA Observations 

The VSA is a CMB interferometer, situated at an altitude 
of 2400 m at Teide Observatory, Tenerife. It has been used 
to measure the CMB angular power spectrum for multipole 
values in the range 150 ^ £ ^ 1500 (Scott et al. 2003; Grainge 
et al. 2003; Dickinson et al. 2004). 

The observations of the dust feature G159.6-18.5 were 
performed during the period September 2005 to July 2008 at 
33 GHz with the "super-extended" array configuration^ . Ta- 
ble 1 provides a summary of the super-extended VSA char- 
acteristics. Observations of calibration sources have shown 
that the VSA interferometric pointing accuracy is better 
than 1'. 

Due to its size 2° x 1°), G159.6-18.5 was observed 
with 11 different pointings, each with a primary beam of 1°.2 
full width at half maximum (FWHM). Each of the 11 point- 
ings were repeated numerous times (Table 2), and were then 
combined to provide the best possible ^^i;- cover age. Table 2 
provides a summary of each of the 11 pointings, including 
the total number of observations, and the corresponding to- 
tal observing time, for each pointing. The rms noise values 
in Table 2 will be discussed in Section 3.2. 

Fig. 2 shows the typical iz, -u-coverage for a 6 hr VSA 
observation with no flagging performed, and one can see 

^ Improved Reprocessing of the IRAS Survey (Miville-Deschenes 
& Lagache 2005). 

The "super-extended" array configuration was the third and 
final array configuration for the VSA. Previously the "compact" 
and "extended" array configurations were in operation. 



the effect of the short spacing problem: the central "hole" 
in the ^^i;- coverage. A lack of short spacings occur because 
there is a finite limit to the smallest aperture spacing (or 
baseline), and this results in leaving a "hole" in the centre 
of the u,v-coverdige. Therefore, the VSA is restricted to only 
observing structure on scales ^ 10 - 40^ and hence our ob- 
servations are "resolving out" any larger scale emission. To 
alleviate this problem and observe this larger scale struc- 
ture, one would need to use single dish observations with 
a ^ 40' beam. 



3 VSA DATA 

3.1 Data Reduction and Calibration 

The VSA data reduction and calibration were performed 
using the software tool reduce. This software was specif- 
ically designed for use with VSA data and is described in 
detail in Dickinson et al. (2004) and references therein. Flux 
density calibration was applied using Tau A with a flux 
density of 350 Jy at 33 GHz. This value was determined 
from observations with the VSA "extended" array (Hafez 
et al. 2008) using Jupiter as the "absolute" calibrator with 
Tb = 146.6 =b 0.75 K at 33 GHz based on WMAP observa- 
tions (Hill et al. 2009). Tau A was chosen due to its close 
proximity (~ 26° away on the sky), which meant that both 
the target and calibrator were available at the same time 
of day. Short observations of the calibrator (?^ 10 mins) al- 
lowed us to use Tau A as both a flux density and phase 
calibrator. Errors due to phase drifts and atmospheric/gain 
errors are typically less than a few percent. Throughout the 
paper we have assumed a conservative absolute calibration 
error of 5 %. 

In addition to performing the flux density and phase cal- 
ibration, REDUCE automatically performs much of the data 
flagging by removing antennas and baselines that were noisy 
or not working correctly. REDUCE also allows the user to fil- 
ter out any contaminating signals such as bright sources 
(e.g. Sun, Moon etc.) passing through the side lobes of the 
VSA beam. This fringe-rate filtering is described in detail 
in Watson et al. (2003). 

After this filtering had been performed, the data were 
fringe-rotated to the centre of the field, and calibrated. All 
the data were inspected and reduced individually, and a u,v 
visibility file for each observation was produced. 

3.2 VSA Map 

In AIPS^ ah the u,v data for each pointing were combined 
to form a single u,v visibility file for each of the 11 point- 
ings. A clean map was then created for each pointing using 
the AIPS task IMAGR, which uses a CLEAN based algo- 
rithm (Hogbom 1974). A small loop gain (0.05) was used in 
the cleaning process to help retain the extended emission, 
in which we are interested. Typical values of the rms noise 
in these images is p:^ 12 - 30 mJy beam~^ (Table 2). These 
rms values were calculated outside the primary beam for 
each image, and show that the images are not dominated 
by thermal noise, but are actually dominated by residual 

^ Astronomical Image Processing System. 
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Table 2. Summary of the VSA observations. 









No. of 


Approx. Total 




Pointing 


HA 


DEC 


observations 


(Jbs. lime 


H.lVl.b. INoise 




(JiUUU) 


(J 2000) 




(hr) 


(mjy beam~^) 


A 


03:41:24 


+32:03:28 


14 


50 


13.3 


B 


03:46:00 


+32:03:28 


8 


20 


16.4 


C 


03:44:12 


+32:03:28 


14 


30 


14.1 


D 


03:42:48 


+32:39:50 


19 


65 


16.4 


E 


03:40:00 


+32:39:50 


9 


25 


19.7 


F 


03:38:36 


+32:03:28 


4 


25 


13.9 


G 


03:40:00 


+31:27:06 


5 


30 


12.0 


H 


03:42:48 


+31:27:06 


7 


25 


14.0 


I 


03:33:32 


+31:34:33 


14 


45 


19.3 


J 


03:37:12 


+32:39:50 


6 


25 


30.6 


K 


03:37:12 


+31:27:06 


5 


20 


26.8 




o 
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Figure 2. Typical w^-u- cover age for a 6 hr VSA observation with 
no flagging performed. Note the central "hole" in the w,t'-coverage 
due to the short spacing problem. X and Y axis are in units of A. 



final VSA map is shown in Fig. 3a with 5 features, all with a 
high signal-to-noise ratio, identified (Al, A2, A3, B and C). 
The central feature, D, in Fig. 3a has a much lower signal- 
to-noise ratio, and as such was not included in our analysis 
but it will be discussed in Section 4.2. 

The IRIS 100 /im image of G159.6-18.5 (Fig. 3b) mostly 
traces the big grain dust emission in the FIR, where the dust 
shell is observed. One can identify that the shell consists of 
a few main features. Comparing Fig. 3a with Fig. 3b, one 
can see that the features observed at 33 GHz are closely 
correlated with the denser dust clumps in G159.6-18.5, but 
not with the ring structure, which does not appear in the 
VSA map. The spatial structure and their implications will 
be discussed in further detail in Section 4. 

The integrated flux density in each of these 6 features, 
shown in Fig. 3a, was calculated by simultaneously fitting 
a gaussian component and a baseline offset using the AIPS 
task JMFIT. The fitting error from jmfit, and the absolute 
calibration error of the VSA (a conservative 5 % - see Sec- 
tion 3.1) were added in quadrature to determine the final 
error. Table 3 lists the position, integrated fiux density and 
deconvolved angular size of each of the 6 features observed 
with the VSA. 



deconvolution effects. An additional contribution to the rms 
noise in these maps arises from the CMB. To access the lev- 
els of contamination from the CMB on the angular scales of 
the VSA, a simulated CMB map, based on a WMAP Syr 
cosmology, was sampled with the VSA sampling distribu- 
tion; this re-sampling process is described in more detail in 
Section 5.1.1. This provides a realistic level of the contribu- 
tion from the CMB in these maps, and we calculated an rms 
of ~ 5 mJy beam~^. This value contributes to the rms noise 
in the maps, but not to a significant level to provde any 
contaminating effects on the 5 features labelled in Fig. 3a. 

A mosaic of all the pointings was then created using 
the AIPS task ltess, which linearly combines each of the 
pointings and applies a primary beam correction in the im- 
age plane. This primary beam correction was applied out 
to ^ 1° (i.e. ^ 2° FWHM) for each individual image. The 



4 ANCILLARY DATA 

Due to its size and range of physical conditions, the Perseus 
molecular complex is a well-studied region, and therefore 
there are a wide variety of ancillary data available, which we 
used to help with our analysis (Table 4). Fig. 4 displays the 
ancillary data^ overlaid with the VSA contours at 33 GHz, 
to help provide a picture of what different emission processes 
are occurring within the region. 



^ Data were downloaded from the Skyview website 

(http://skyview.gsfc.nasa.gov), the LAMBDA website 

(http://lambda.gsfc.nasa.gov) and the COMPLETE Survey 
website (http://www.cfa.harvard.edu/COMPLETE). 
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Figure 3. G159.6-18.5 as observed with the VSA at 33 GHz (left) and IRIS 100 fim (right). The location of the dust sheh (dashed line), 
the O9.5-B0V star, HD 278942 (cross), and the position of the 6 features in the VSA image (Al, A2, A3, B, C and D) are indicated. The 
5 features in the VSA map, with high signal-to-noise ratios (Al, A2, A3, B and C), are discussed and analysed throughout the rest of the 
present paper. The central feature, D, with a much lower signal-to-noise ratio, was not included in our analysis, but will be discussed in 
Section 4.2. The VSA ^ 7' FWHM synthesized beam is shown in the bottom left-hand corner of the image, and the contours correspond 
to 10, 20, 30, 40, 50, 60, 70, 80, 90 % of the peak flux which is 200 mJy beam"!. The units of the IRIS image are MJy sr-^. 



Table 3. Position, integrated flux density and deconvolved angular size of the 6 features observed with the VSA at 33 GHz. Features 
Al, A2, A3, B and C are analysed throughout the rest of the present paper, but feature D was not included in our analysis due to the 
much lower signal-to-noise ratio. 



Feature 


RA 


DEC 


Flux Density 


Major Axis 


Minor Axis 




(J2000) 


(J2000) 


(Jy) 


(arcmin) 


(arcmin) 


Al 


03:44:30 


+32:11:27 


0.53 ± 0.12 


21.6 ± 4.4 


4.4 ± 3.8 


A2 


03:43:27 


+32:00:17 


1.39 ± 0.20 


24.3 ± 3.1 


12.0 ± 2.2 


A3 


03:43:14 


+31:46:25 


0.37 ± 0.13 


16.1 ± 5.8 


9.8 ± 5.4 


B 


03:40:41 


+31:16:31 


1.23 ± 0.21 


20.5 ± 4.8 


17.6 ± 4.8 


C 


03:37:49 


+31:27:57 


0.86 ± 0.14 


14.0 ± 3.1 


13.5 ± 2.9 


D 


03:40:49 


+31:54:38 


0.18 ± 0.10 


16.0 ± 9.1 


4.9 ± 5.7 



4.1 Radio Frequency Data 

At frequencies below ~ 60 GHz, the emission is expected 
to be dominated by synchrotron and free-free emission. 
However, since G159.6-18.5 is far from the Galactic plane 
(/ = 159°. 6 and h — —18° 5) and is not in the vicinity of 
any known supernova remnants, the synchrotron emission 
is expected to be negligible at 33 GHz. Therefore the low 
frequency ancillary data help us to understand the free-free 
emission. 

In the Stockert 1.4 GHz survey (Fig. 4a), one can iden- 
tify a region of diffuse emission within the shell that is offset 
from the centre by ~ 10 - 15^ This diffuse region of ~ 40 ar- 
cmin diameter has also been mapped at 2.7 GHz by Reich 
& Reich (2009). There appears to be little or no emission at 
1.4 GHz from the 5 features observed at 33 GHz with the 
VSA. This agrees with the current understanding of the re- 
gion, which is believed to contain a bubble of Hii gas within 
the dust shell (see Section 2.1). The Ha image (Fig. 4b) ap- 
pears to confirm this hypothesis, with the Ha emission also 
being confined to the interior of the shell, and also offset 



similar to the Stockert observations, which we would expect 
as it traces the hot (Te ~ 10^ K) ionized gas. 

The NRAO VLA Sky Survey (NVSS) at 1.4 GHz and 
the GB6 4.85 GHz survey (Fig. 4c and Fig. 4d respectively) 
have much higher resolution than both the Stockert survey 
and the VSA observations, and are very useful for identify- 
ing compact radio sources, including Hii features. We find 
that there are three radio sources in the G159.6-18.5 re- 
gion: 4C 31.14 to the south-east of the ring; 3C 92 within 
the ring; and 4C 32.14 to the north-west of the ring. Each 
of these sources, 4C 31.14, 3C 92 and 4C 32.14, are extra- 
galactic with spectral indices^ of a ~ —0.7, —1.0 and —0.03 
respectively. The flat spectrum quasar 4C 32.14 was removed 
from the VSA map by subtracting a CLEAN component 
model from the UV data. All 5 features observed with the 
VSA, however, appear to be free from any source contami- 
nation. Extra-galactic source counts at 15 GHz (Waldram et 



^ Throughout this paper, the flux density spectral index is de- 
flned as Sjy oc u^. 
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Table 4. Summary of the ancillary data used in our analy- 
sis. References are: [1] Reich & Reich (1986); [2] Finkbeiner 
(2003); [3]: Condon et al. (1998); [4] Condon, Broderick & Seiel- 
stad (1989); [5] Ridge et al. (2006a); [6] Hinshaw et al. (2009); 
[7] Miville-Deschenes & Lagache (2005). 



Telescope/ 


Frequency/ 


Angular 


Ref. 


Survey 


Wavelength 


Resolution 




Stockert 


1.4 GHz 


34' 


[1] 


WHAM/VTSS/SHASSA 


Ha 


6' 


[2] 


NVSS 


1.4 GHz 


0'.75 


[3] 


GB6 


4.85 GHz 


3^.5 


[4] 


2MASS/NICER 


Av 


5' 


[5] 


WMAP 


94 GHz 


12^.6 


[6] 


IRAS/IRIS 


12 /im 


3^.8 


[7] 


IRAS/IRIS 


25 /im 


3^.8 


[7] 


IRAS/IRIS 


60 /im 


4^0 


[71 


IRAS/IRIS 


100 /xm 


4^3 


[7] 



al. 2003) imply that the 30 GHz extra-galactic source popu- 
lation cannot remotely account for the 5 observed features. 
Also the fact that the features we observe are extended, 
strongly rules out the possibility of significant extra-galactic 
source contamination. 

Weak large scale emission, predominantly resolved out, 
can be seen within the shell in the GB6 image, and it is 
aligned with the Stockert survey observations. However, this 
large scale diffuse emission does not extend as far as the 5 
features observed with the VSA. 

The visual extinction (Av) map (Fig. 4e) shows that the 
extinction is quite high throughout the entire region 2 - 
10 mag), suggesting a large quantity of dust in the vicinity 
(comparing Fig. 4e with Fig. 4f confirms that the extinction 
is due to the large, cold dust grains with T < 15 K, which are 
traced by the 94 GHz emission) . This large quantity of dust 
makes it very difficult to use the Ha emission as a means 
to calculate an estimate of the free-free emission within the 
region (Dickinson et al. 2003). The reason for this is that it is 
very difficult to judge what fraction of the dust, causing the 
extinction, is situated between us and the molecular cloud. 

However, to establish a rough estimate of whether the 
emission in the VSA map could be due entirely to free-free 
emission or not, we calculated the associated Ha emission 
you would expect, assuming that all the emission observed 
with the VSA at 33 GHz was due to free-free emission, based 
on the assumptions of Dickinson et al. (2003). Using the 
conservative assumption that all the extinguishing dust is 
located between us and G159.6-18.5, an upper limit esti- 
mate of the actual Ho; emission was calculated. We found 
that the Ha emission is much too faint (typically on average 
by a factor of 20) to account for the observed 33 GHz emis- 
sion. This supports the idea that the emission observed at 
33 GHz is not originating entirely from free-free emission. 
In Section 5.1.1 we present a more rigorous analysis that 
constrains the level of free-free emission observed with the 
VSA. 

4.2 Higher Frequency Data 

Emission at higher frequencies (^60 GHz) is dominated by 
thermal dust emission, and we used the WMAP 5-yr total 



Table 5. Upper limits to the free-free emission in the 5 features. 
A 4.85 GHz 3a flux density upper limit based on the GB6 image. 
This flux density is scaled to 33 GHz using a nominal free-free 
spectral index of a = —0.12. The last column represents the per- 
centage of free-free emission observed with the VSA at 33 GHz. 



Feature 


GB6 

S4.85 GHz 
(mJy) 


GB6 

S33 GHz 
(mJy) 


Filtered VSA 

S33 GHz 
(mJy) 


% due to 
free-free 
(3a) 


Al 


< 68.3 


< 54.3 


178 lb 52 


< 31 


A2 


< 38.8 


< 30.8 


155 di 41 


< 20 


A3 


< 32.9 


< 26.2 


167 lb 56 


< 16 


B 


< 41.0 


< 32.6 


354 ± 69 


< 9 


G 


< 43.1 


< 34.2 


95 ± 38 


< 36 



intensity W-band map, and the IRIS 100, 60, 25 and 12 /xm 
maps to enable us to understand this emission. 

The WMAP W-band (94 GHz) image (Fig. 4f) traces 
the large, cold dust grains (T < 15 K) and shows that there 
is emission on large angular scales in the central region of 
G159.6-18.5. 

The IRIS images in Fig. 5 display the emission in the 
mid-IR to FIR, and highlight the high degree of correlation 
between the IR and the microwave emission. 

The region in the centre of the shell is interesting as 
there appear to be two distinct FIR components, offset 
by ~ 20^ One component to the east, which is visible in 
all 4 IRIS bands, and another more central component that 
appears to peak at ^ 60 /im, suggestive that this is much 
warmer (T ^ 70 K) than the surrounding dust (T p:^ 30 - 
40 K as displayed on the plots in Fig. 6). From Fig. 5c, one 
can see that the location of the 09.5-BOV star, HD 278942, 
corresponds to the brighter central component, and hence 
this could be acting as a heat source for the dust, and there- 
fore producing the higher dust temperatures observed lo- 
cally. 

These IRIS images and the correlation between the mi- 
crowave and the FIR emission will be discussed further in 
Section 5.1.2 and Section 6.3 respectively. 

5 SPECTRAL ENERGY DISTRIBUTIONS OF 
THE 5 REGIONS 

To identify whether the emission observed in the VSA map 
shown in Fig. 3a is actually anomalous, we looked at the 
spectrum to determine whether there was an excess over 
the level of free-free and thermal dust emission. In order to 
constrain the levels of free-free emission we used the GB6 
4.85 GHz observations (Section 5.1.1), and to constrain the 
levels of thermal dust emission we used the IRIS observa- 
tions in the FIR regime (Section 5.1.2) and the WMAP W- 
band data in the Raleigh-Jeans regime (Section 5.1.3). 

5.1 Integrated Flux Density Spectra of Individual 

Features 

5.1.1 GB6 

As discussed in Section 2.2, the VSA only observes struc- 
tures on angular scales of ?^ 10 - 40^ Therefore, to deter- 
mine the level of free-free emission in the VSA map, we need 
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Figure 4. Ancillary data overlaid with the VSA contours and synthesized beam, as shown in Fig. 3a, displaying the different emission 
mechanisms in grey scale occurring in the region: a) Stockert 1.4 GHz in units of mK (T^); b) Ho; in units of R; c) NVSS 1.4 GHz in 
units of mJy beam"-*^; d) GB6 4.85 GHz in units of mJy beam"-*^; e) Av in units of mag; f) WMAP 94 GHz in units of Jy beam"-*^. The 
dashed line displays the position of the dust shell, and the cross marks the position of the central star, HD 278942. 
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Figure 5. False colour IRIS images of G159.6-18.5 processed with the VSA sampling distribution (see Section 5.1.2) overlaid with the 
VSA contours and synthesized beam, as shown in Fig. 3a: a) IRIS 100 jim in units of kJy beam"-*^; b) 60 jim in units of kjy beam"-*^; 
c) 25 fim in units of Jy beam"-*^; d) 12 /xm in units of Jy beam"-*^. The dashed line displays the position of the dust shell, and the cross 
marks the position of the central star, HD 278942. These images highlight the correlation between the emission at microwave and IR 
wavelengths. 



some low frequency observations with a similar angular reso- 
lution, namely, observations with an equivalent range of u,v- 
coverage. In Section 4.1 we discussed the three available low 
frequency surveys: Stockert (34' FWHM at 1.4 GHz); NVSS 
(0'.75 FWHM at 1.4 GHz); GB6 (3'.5 FWHM at 4.85 GHz). 
None of these surveys provide an ideal match for the VSA. 
We cannot use the Stockert survey due to the lack of short- 
spacings in the VSA coverage, where all the large scale struc- 
ture has been resolved out and the NVSS also suffers from 
an absence of short-spacings, which means that we cannot 
use it. This leaves the GB6 survey. 

The original GB6 observations have been filtered to re- 
move all structure on scales larger than ^ 20' (Condon, 
Broderick & Seielstad 1989) and hence cannot be used di- 
rectly to assess the levels of free-free emission in the VSA 
image. To allow us to make use of the GB6 image, we had to 
filter the VSA image so as to harmonise the u,v ranges sam- 



pled. This involved removing all the short (^ 172 A) VSA 
baselines. The VSA sampling distribution was then applied 
to the GB6 image using the AIPS task UVSUB. This pro- 
duces a u,v visibility file that has been sampled with the 
VSA sampling distribution, but with the GB6 visibilities. 
This was then CLEANed using the AIPS task IMAGR in a 
similar manner to the VSA image, creating a new GB6 im- 
age. The filtered VSA image was then directly comparable 
to this newly re-sampled GB6 image, and hence allowed us 
to determine the contribution of free-free emission on these 
angular scales. 

The flux density of the 5 features in the filtered VSA 
image was calculated using jmfit applying both a gaussian 
and baseline component and is given in Table 5). As dis- 
cussed in Section 4.1, all the features in the VSA image 
appear to show no detectable emission in the GB6 image, 
therefore, using the results of jmfit on the filtered VSA 
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Table 6. Integrated flux density for each of the 5 features identifled in the VSA map at 4.85 GHz, 33 GHz, 94 GHz, 100 /xm and 60 jim. 
All errors quoted include a fitting error and an absolute calibration error (5 % for VSA, 5 % for WMAP data and 10 % for the IRIS 
data) combined in quadrature. "''Flux density at 4.85 GHz is a conservative 3a upper limit based on the GB6 image (see Section 5.1.1 
for details). 



Feature 


S4.85 GHz^ 


S33 GHz 


S94 GHz 


SlOO fxm 


SeO iJ,m 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


Al 


< 0.20 


0.53 ± 0.12 


0.65 ± 0.63 


5000 ± 560 


3200 ± 330 


A2 


< 0.35 


1.39 ± 0.20 


1.99 ± 0.89 


2200 ± 460 


810 ± 160 


A3 


< 0.07 


0.37 lb 0.13 


0.34 di 0.55 


770 lb 260 


260 lb 90 


B 


< 0.14 


1.23 lb 0.21 


0.75 lb 1.48 


2500 lb 520 


850 lb 170 


C 


< 0.39 


0.86 lb 0.14 


< 1.49 


1800 lb 440 


330 lb 110 



image, along with the rms noise in the same region of the 
re-sampled GB6 image, a 3a upper limit for the integrated 
flux density at 4.85 GHz, in each feature, was determined. 
This is a conservative upper limit since we are dealing with 
integrated fluxes rather than with peak surface brightnesses. 

Table 5 lists the integrated flux density for the 5 fea- 
tures at 4.85 GHz and 33 GHz, calculated as described pre- 
viously in this section. The 3a flux density upper limit at 
4.85 GHz has been scaled to 33 GHz, using a nominal free- 
free spectral index of o; = —0.12 (Dickinson et al. 2003). A 
direct comparison between the emission found in the GB6 
and VSA images is given as a percentage, in the last column 
in Table 5. This suggests that the contribution of free-free 
emission at 33 GHz is far from dominant, and rules out the 
possibility that the emission in the VSA map is due entirely 
to free-free emission. By assuming that this ratio between 
the 4.85 GHz and the 33 GHz emission is independent of 
angular scale, this allows us to compute an estimate for the 
contribution of free-free emission at 33 GHz on the angular 
scales of the unfiltered VSA image (see Table 6). 

5.1.2 IRIS 

The 4 IRIS images at 100, 60, 25 and 12 /xm have compa- 
rable resolution to the VSA (Table 4). The VSA sampling 
distribution was applied to the IRIS images using the AIPS 
task UVSUB (in a similar manner to the GB6 images as de- 
scribed in Section 5.1.1). These new IRIS images (see Fig. 5) 
can then be directly compared with the VSA image, jmfit 
was performed on these images to calculate the integrated 
flux density in each of the 5 features in the VSA map. This 
was performed using the same JMFIT parameters obtained 
from the VSA image, and the fitting error and the absolute 
calibration error (conservatively assumed to be « 10 %) were 
combined in quadrature to determine the overall error (see 
Table 6). 

5.1.3 WMAP W-hand 

The 5-yr WMAP W-band map was used to constrain the 
Raleigh- Jeans tail of the thermal dust emission. As dis- 
cussed in Section 4.2 most of the emission at 94 GHz oc- 
curs within the central region of G 159. 6-1 8. 5, and extends 
to IC 348. To calculate the integrated fiux density of each of 
the 5 features at 94 GHz, the VSA image was smoothed to 
a 12'.6 resolution, and then using JMFIT each of the features 
was fitted with a gaussian component. The size and position 



of these gaussians were then kept fixed, and fitted for in the 
WMAP W-band image. This accounted for the difference in 
angular resolution between the two images. 



5.2 SED Determination 

Fig. 6 shows the fiux density spectrum for each of the 5 
features identified in Fig. 3a. In Section 4 we discussed the 
emission mechanisms expected to occur within this region 
using the ancillary data available. Free-free emission occurs 
at the radio frequencies, and thermal (vibrational) dust is 
expected to occur at FIR wavelengths. 

We determined an upper limit for the free-free emis- 
sion at 4.85 GHz, and since we know that a typical free- 
free emission spectrum follows a power-law with spectral in- 
dex a = —0.12, we plotted this normalized to our upper limit 
at 4.85 GHz. We note that at frequencies below ~ 1 GHz, the 
free-free emission can become optically thick (r > 1) with a 
spectral index oia^ -\-2^ but at higher frequencies, the free- 
free emission becomes optically thin (r <^ 1) and follows the 
standard power-law with spectral index a = —0.12. 

For the thermal dust emission, we used a modified 
black-body curve function, z/^'^^^*^^B(z/,Tdtist), where jSdust 
is the dust emissivity index and B(i/,Tdnst) is the black-body 
function for a dust temperature, Tdust- A single tempera- 
ture dust component, modelled by this modified black-body 
curve, was fitted to the 94 GHz, 100 /xm and 60 /xm data 
points for each feature individually. The best fit values for 
both j3dust and Tdust are displayed on each of the plots in 
Fig. 6. 

Also plotted in Fig. 6 are the 25 and 12 /xm data points. 
These were not fitted by the single temperature modified 
black-body curve because the emission at these wavelengths 
is not due to the same dust grain population as the 100 and 
60 fim.. The 100 and 60 fim. emission is dominated by BGs, 
while the 25 and 12 /xm emission is dominated by VSGs. The 
12 /im emission has an additional contribution from PAHs. 



6 DISCUSSION 
6.1 SED Modelling 

In each plot in Fig. 6, it is possible to identify an excess of 
emission present at 33 GHz that cannot be explained by ei- 
ther free-free or thermal dust emission. Table 7 displays this 
fraction of excess emission at 33 GHz, and the corresponding 
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Figure 6. Integrated flux density spectra for the 5 features shown in Fig. 3a. A power-law with spectral index a = —0.12 is plotted to 
represent the free-free emission 3<t upper limit (dotted line) and a modifled black-body curve, fitted to the 94 GHz, 100 yLtm and 60 /Ltm 
data points, is plotted to represent the thermal dust emission (dashed line). The best fitting values for Tdust and l3dust are displayed on 
each of the plots. For feature C, the 94 GHz data point was only an upper limit, and therefore the modified black-body curve is fitted 
with a fixed dust emissivity value of jS^ust — 1-5. The VSA 33 GHz data point is not accounted for by either of these curves and is fitted 
by a linear combination (0.5CNM + 0.5MC) of Draine & Lazarian (1998) spinning dust models (dash-dot line). Also plotted are the 
25 /Ltm and 12 yLtm data points, but these were not included in the SED fitting procedure. 
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Table 7. Fraction of excess emission at 33 GHz, and the signifi- 
cance, in each of the five features. 



Feature Fraction of Significance 
excess emission 

Al > 0.69 > 3.0(7 

A2 > 0.79 > 5.6a 

A3 > 0.82 > 2.5(1 

B > 0.90 > 5.2(7 

C > 0.62 > 3.8(7 



significance, for each of the five features, confirming that a 
substantial fraction of the emission at 33 GHz is anomalous. 

For features Al, A2, A3 and B, the 94 GHz data point 
provides a constraint on the shape of the thermal dust 
curve, providing the best fitting values for the dust tem- 
perature Tdust, and the dust emissivity index, f3dust (see 
plots in Fig. 6). However, for feature C, the 94 GHz data 
point is an upper limit and hence does not provide a reli- 
able constraint on the thermal dust curve. Therefore, for this 
feature, we used a fixed value of (3dust = 1-5. This is a typ- 
ical value of pdust (Dupac et al. 2003), and to ensure that 
selecting this value was reasonable, we investigated using 
a value of f3dust = 1-1 (towards the lowest measured value 
for Pdust (Dupac et al. 2003)), and found that this only re- 
duced the fraction of excess emission at 33 GHz from 62 % 
to 48 %. This confirms that there is an excess of emission 
observed with the VSA at 33 GHz that cannot be explained 
by either free-free or thermal dust emission. Therefore to fit 
the 33 GHz data point, a peaked spectrum is required. 

Peaked spectra can be caused by emission from ultra- 
compact (UC) Hii regions, gigahertz-peaked spectra (GPS), 
magnetic dipole emission or electric dipole emission. UCHll 
regions are very compact, optically thick, Hii regions while 
GPS are radio sources at high redshifts with synchrotron 
self-absorption. A detailed search of the IRAS point source 
catalogue (Beichman et al. 1988) revealed that there are only 
a few sources (< 4 %), within 1° of the centre of our image, 
that meet the colour criterion set by Wood & Churchwell 
(1989a) for UCHil regions. Further analysis reveals that to 
fit the data with emission from an UCHii region would re- 
quire an emission measure of ^ 2 x 10^ cm~^pc and an 
angular size of ^ 1'' . These parameters are towards the up- 
per end of the observed values for UCHii regions (Wood 
& Churchwell 1989b), and also imply that UCHll regions 
would appear as point-like sources in our image, which is 
not the case. We would also expect GPS to appear as point- 
like sources in our image, and to be visible in the NVSS and 
GB6 surveys discussed in Section 4.1. Neither of these are 
true. 

Magnetic dipole emission has been ruled out due to 
polarization observations of the region by Battistelli et al. 
(2006). 

Therefore, having ruled out UCHii regions, GPS and 
magnetic dipole emission as viable options, we have fit- 
ted a Draine & Lazarian (1998) spinning dust model for 
electric dipole radiation. Since we only have the one data 
point at 33 GHz, we fitted a linear combination of the 
models for the cold neutral medium and the dense molecu- 
lar cloud (0.5CNM + 0.5MC) to ensure the peak occurred 



Table 8. Dust emissivities at ^ 30 GHz relative to 100 /xm for the 
5 features observed in this work, and 6 other sources to provide a 
comparison. The 6 Hn regions is a mean value based on 6 southern 
Hii regions, LPH 96 is also an Hii region and LDN 1622 is a dark 
cloud. The all-sky value is a mean value outside the Kp2 mask, 
and the 15 regions value is a mean value of 15 high latitude regions 
both from WMAP. G159.6-18.5 is the value obtained with the 
COSMOSOMAS observations. 



Source 


Dust Emissivity 
ijiK (MJy sr-l)-l 


Reference 


Al 


2.8 lb 0.7 


This paper 


A2 


16.4 lb 4.1 


This paper 


A3 


12.8 lb 6.1 


This paper 


B 


13.2 ± 3.6 


This paper 


C 


13.0 ± 3.2 


This paper 


6 Hn regions 


3.3 ± 1.7 


Dickinson et al. (2007) 


LPH96 


5.8 ± 2.3 


Dickinson et al. (2006) 


All-sky 


10.9 ± 1.1 


Davies et al. (2006) 


15 regions 


11.2 ± 1.5 


Davies et al. (2006) 


LDN1622 


24.1 ± 0.7 


Casassus et al. (2006) 


G159.6-18.5 


15.7 lb 0.3 


Watson et al. (2005) 



at 33 GHz. Clearly more data in the frequency range 5 - 
60 GHz are required to confirm this fit. 

6.2 Dust Emissivities 

The dust emissivity, the ratio between the microwave and 
IR emission, provides a convenient way^ to represent the 
excess emission found in Fig. 6 and is measured in units 
of fiK (MJy sr~^)~^. Using the values in Table 6, the dust 
emissivity between 33 GHz and 100 /im was calculated for 
each of the 5 features and are tabulated in Table 8. Also 
listed in Table 8 is the dust emissivity for 6 other sources: 
2 values based on Hll regions and 4 values for typical high 
Galactic latitudes. 

The radio emissivity of dust at ^ 30 GHz is known 
to have a typical value of ^ 10 /jK (MJy sr~^)~^ with a 
factor of 2 variation at high Galactic latitudes (Davies et 
al. 2006). The dust emissivity for the 5 features analysed in 
this work are consistent with this value. Comparing the dust 
emissivities listed in Table 8, we find that the features A2, 
A3, B and C have a dust emissivity similar to the typical 
value found for high Galactic latitudes, while feature Al, 
with a much lower dust emissivity, is similar to that for Hll 
regions. 

The lower dust emissivity per unit of 100 /xm surface 
brightness in Hll regions (see Table 8) may be the result of 
the relative difference in physical environmental conditions 
compared with those in the general ISM. For example, the 
warmer dust resulting from the strong UV radiation field in 
an Hii region will emit more strongly at 100 /xm, thereby re- 
ducing the radio to FIR emissivity ratio for a given amount 
of dust. Also the more energetic radiation field may disas- 
sociate the VSGs/PAHs responsible for the spinning dust 
emission. Feature Al, associated with IC 348, has an emis- 
sivity similar to an Hii region although no Hii region with a 

^ The dust emissivity is calculated relative to the 100 yLtm emis- 
sion, and is therefore dependent on the temperature of the dust. 
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Figure 7. Pearson correlation coefficient plotted as a function of 
wavelength. A fit to all 4 IRIS wavebands (dotted line), and a 
fit to the 25, 60 and 100 /xm wavebands (dash-dot-dot line) are 
plotted. Both fits show a significant non zero gradient, at a 2.5a 
and 3.8(7 level respectively, with the correlation increasing as the 
wavelength decreases, la errorbars are plotted. 



flux density greater than 5 mJy at radio frequencies is seen 
with the VSA at the position of IC348. Nevertheless, there 
is a significant UV radiation field arising from the low mass 
stars comprising the open cluster in the IC348 reflection neb- 
ulosity and there may additionally be some UV illumination 
from the nearby Perseus 0B2 association. The overall radi- 
ation field in IC 348 is 10 - 100 times that of the general 
ISM (Bachiller et al. 1987). The higher dust temperature 
{T^dust ~ 40 K) compared with that {Tdust ~ 30 K) in the 
cooler features confirms the presence of this higher radiation 
field in the Al region. The lack of a detectable Hll region 
indicates a less intense radiation field or a lower gas density 
than in a normal Hli region. 

Ignoring feature Al, with its lower dust emissivity, 
the mean dust emissivity for features A2, A3, B and 
C (13.8 =b 2.2) is consistent with the integrated value cal- 
culated for the whole region by Watson et al. (2005). These 
comparisons of dust emissivity confirm that our detection of 
excess emission is in agreement with previous observations, 
and can be explained in terms of spinning dust. 

6.3 Correlation between the Microwave and IR 
Observations 

A visual comparison of the IRIS images overlaid with the 
VSA contours in Fig. 5 suggests a possible correlation be- 
tween the microwave and IR emission. A Pearson correla- 
tion analysis was performed, and Fig. 7 illustrates how good 
the overall correlation is 0.8), and displays the calcu- 
lated Pearson correlation coefficient as a function of IRIS 
wavelength. The plotted error bars are la, and were deter- 
mined using the Fisher z transformation, which converts 
the non-normal Pearson sampling distribution to a normal 
distribution with standard errors. 

Also plotted on Fig. 7 are 2 fits: one fit applied to all 4 
IRIS wavelengths (12, 25, 60 and 100 /xm), and one fit ap- 
plied to only 3 of the IRIS wavelengths (25, 60 and 100 /im). 
Both fits display a non-zero gradient at a 2.5cr level and 
a 3.8cr level respectively. This implies that the correlation 



between the microwave and the IR increases as the wave- 
length decreases. The shorter IR wavelengths are known to 
be due to emission from the smallest dust grains (VSGs and 
PAHs), and hence the fact that the correlation is stronger 
with shorter wavelengths reinforces the hypothesis of this 
anomalous emission being caused by spinning dust. 

This correlation analysis was performed across 4 of the 
5 features in the VSA image (A2, A3, B and C). Feature 
Al was excluded because it has different physical conditions 
compared to the other 4 features as shown in Section 6.2 
and Table 8. When feature Al is excluded from the corre- 
lation analysis, the overall correlation coefficient increases 
by ^ 4 %. 

The 12 fim data point in Fig. 7 appears to be low, 
however we must remember that although the IRIS data 
are an improved version of the original IRAS data, there 
is still some substantial "contamination" from point sources 
in these images. These point sources would be expected to 
dominate at the shorter wavelengths and hence this could 
explain the apparent decrease in correlation at 12 /xm. 

Further investigation using data from the Spitzer Space 
Telescope will be performed in a follow-up paper. This 
should clear up our understanding of whether the point 
source contamination is responsible for the apparent peak 
at 25 fim. 



6.4 Anomalous Emission and Angular Scale 

G159.6-18.5 was observed by Watson et al. (2005) with the 
COSMOSOMAS experiment on angular scales of ^ 1^, and 
hence could not resolve any of the 5 individual features ob- 
served with the VSA in this present work. Using both results 
allow us to investigate how the anomalous emission varies 
with angular scale. The total WMAP flux density measured 
by Watson et al. (2005) was 40.3 =b 0.4 Jy at 33 GHz, and 
the total flux density in the 5 features in the VSA map 
was found to be 4.4 zb 0.4 Jy (Table 6). This implies that 
the VSA is only observing a small fraction (?^ 10 %) of the 
anomalous emission observed by COSMOSOMAS, and that 
the bulk of the anomalous emission is originating from some 
large scale structure that is being resolved out by the VSA, 
due to the short-spacing problem. This idea also explains 
why the integrated flux density in the flltered VSA image, 
0.96 =b 0.12 Jy (Table 5) is much lower than the flux den- 
sity in the unfiltered VSA image. Therefore, the features 
observed by the VSA must be "peaks" of anomalous emis- 
sion, located on a large scale, diffuse cloud responsible for 
the 33 GHz emission inWMAP. 



7 CONCLUSIONS 

Observations of G159.6-18.5 with the VSA in its super- 
extended array configuration with a synthesized beam 
of ~ 7' FWHM provide clear evidence for excess emission 
at 33 GHz in 5 distinct dust-correlated components. This 
anomalous emission cannot be explained by either free-free 
or thermal dust emission, and as such produces a peaked 
spectrum. After ruling out both UCHll regions and GPS by 
observations of ancillary data, and magnetic dipole emission 
based on previous observations (Battistelli et al. 2006), we 
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explained the excess emission by fitting for electric dipole ra- 
diation (spinning dust) as predicted by Draine & Lazarian 
(1998). 

From the fitted SEDs, feature Al was found to 
have Tdust ~ 40 K, which was significantly hotter than the 
other 4 features with Tdust ~ 30 K. The central component 
of the shell peaks at ^ 60 //m, and as such was found to 
have Tdust ^ 70 K. We suspect this central heating is due 
to the presence of the central 09.5-BOV star, HD 278942. 

The dust emissivity for each of the 5 features were cal- 
culated, and found to be in agreement with the values from 
previous observations where anomalous emission has been 
detected. Feature Al was found to have a similar value to 
that of Hii regions, while the other 4 features all agree with 
the typical value for high Galactic latitudes. 

By performing a correlation analysis between the VSA 
and IRIS observations, we have shown that the correlation is 
stronger at shorter wavelengths, suggesting that this anoma- 
lous emission is caused by the smallest dust grains such as 
VSGs or PAHs. It is these dust grains, which current spin- 
ning dust models (Draine & Lazarian 1998; Ali-Haimoud et 
al. 2009) predict as the source of the electric-dipole radia- 
tion, and hence these correlation results imply that what we 
are observing is emission from spinning dust grains. Spec- 
troscopic observations (Iglesias-Groth et al. 2008) have pro- 
vided evidence for the existence of the naphthalene cation (a 
simple PAH) in the G 159. 6-1 8. 5 region, which could be re- 
sponsible for the anomalous emission. 

To gain a better understanding of the origin and nature 
of this anomalous emission, we compared our results with 
the previous analysis of this region performed by Watson 
et al. (2005). This comparison revealed that the bulk of this 
anomalous emission (^ 90 %) appears to be originating from 
some large scale structure, which is being resolved out by the 
VSA i.e. that the anomalous emission in G159.6-18.5 is very 
diffuse and not concentrated in the 5 features observed with 
the VSA. 

Future high-resolution observations at microwave wave- 
lengths, combined with the Spitzer IR data, is needed to 
investigate in more detail the nature of the anomalous emis- 
sion. The Planck and Herschel satellites will provide data in 
the microwave, sub- mm and FIR regimes which are required 
to measure the spectrum of anomalous dust in the frequency 
range 10 - 100 GHz. 
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